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Measuring the local quantum capacitance of graphene using a strongly coupled 

graphene nanoribbon 
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We present electrical transport measnrements of a van-der-Waals heterostructure consisting of a 
graphene nanoribbon separated by a thin boron nitride layer from a micron-sized graphene sheet. 

The interplay between the two layers is discussed in terms of screening or, alternatively, quantum 
capacitance. The ribbon can be tuned into the transport gap by applying gate voltages. Multiple 
sites of localized charge leading to Coulomb blockade are observed in agreement with previous 
experiments. Due to the strong capacitive coupling between the ribbon and the graphene top layer 
sheet, the evolution of the Coulomb blockade peaks in gate voltages can be used to obtain the 
local density of states and therefore the quantum capacitance of the graphene top layer. Spatially 
varying density and doping are found which are attributed to a spatial variation of the dielectric 
due to fabrication imperfections. 

PACS numbers: 71.15.Mb, 81.05.ue, 72.80.Vp 


One of the advantages of layered materials such as 
graphene, hexagonal boron nitride (hBN) or various 
dichalcogenides, is the possibility to fabricate stacks in 
order to obtain a heterostructure material with new prop¬ 
erties. This was recently exploited to fabricate a variety 
of different device structures m and to investigate a 
number of novel effects: localization in super clean sys¬ 
tems laiin], current drag [iiHia and the emergence of 
superlattices [I4lll6j . Most of these experiments were 
performed with so-called “double-layer graphene”, con¬ 
sisting of two graphene layers separated by a thin insula¬ 
tor layer. In contrast to Bernal-stacked bilayer graphene, 
charge transfer between the two layers is not possible. 

Above experiments were conducted with micron-sized 
graphene sheets. In this paper we develop this idea 
one step further: we pattern the lower graphene sheet 
laterally into a nanoribbon (e.g. Refs. UMI]) and 
stack an insulator plus a micron-sized graphene sheet on 
top. This experiment is a crucial step towards fabricat¬ 
ing and understanding stacked graphene nanostructures. 
If the nanoribbon is tuned into the so-called transport 
gap [171 [E] ) Coulomb blockade can be observed [171 HI] ■ 
The localized charge in the ribbon serves as a sensitive 
and local detector for its electrostatic surrounding and 
allows us to probe properties of the graphene top layer. 
This is conceptually similar to an SET above a 2DEG [22] 
and scanning SET measurements [231 HI] • Despite using 
a comparably thick insulator layer (13nm), we can still 
reach the regime where the spacing between the layers 
is smaller than the average distance between two charge 
carriers within one layer [la¬ 
in this paper we present simultaneous electrical trans¬ 
port measurements through a graphene nanoribbon and 
a graphene sheet (top layer) that are separated by a layer 
of hBN, and therefore strongly capacitively coupled. By 
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applying gate voltages, we can independently tune both 
the graphene top layer and the graphene ribbon to their 
respective Dirac point. If the graphene ribbon is tuned 
close to the Dirac point, electrical transport is dominated 
by Coulomb blockade due to localized charges. We use 
the evolution of Coulomb resonances in gate voltages to 
extract the local density and quantum capacitance of the 
graphene top layer. 

The investigated device is shown in Fig. [^ The single 
layer graphene nanoribbon was patterned by reactive ion 
etching, is 200 nm long and about 65 nm wide. The rib¬ 
bon was fabricated on silicon dioxide since using hexago¬ 
nal boron nitride (hBN) substrates complicates the fab¬ 
rication while not significantly changing the electronic 
transport properties for such ribbons |21j . Two graphene 
side gates are located left and right of the ribbon. As vis¬ 
ible in Fig. [gi, there are PMMA residues on the ribbon 
that remained after chemical cleaning. Hereafter, the 
layer with the ribbon and the side gates is called “rib¬ 
bon layer”. On top of the ribbon layer, a hBN flake with 
a thickness of 13 nm - as determined by scanning force 
microscopy (SFM) - was deposited. An additional single 
layer graphene flake was transferred on top of the hBN 
flake (top layer). The top layer flake exhibits wrinkles 
and bubbles and is larger than the ribbon layer graphene 
flake as shown in Fig. [^. As shown in Fig. [^, the area 
directly above the ribbon is free of big wrinkles and bub¬ 
bles. For fabrication details see Refs. [mus]. 

All measurements were performed at a temperature of 
1.3 K and all DC voltages were applied relative to the 
ribbon as shown in Fig. No leakage currents between 
the ribbon, the top layer or any of the gates could be 
measured. Unless stated differently, Vsg,i = Use ,2 = 
OV. 

Fig. [^ shows the resistance of the graphene top layer 
as a function of Vtl and Vbg at a fixed current bias of 
200 nA. The voltage drop is measured over the two con¬ 
tacts marked in Fig. [^ while all other top layer contacts 
are left floating. Contact and cable resistances of a few 
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Figure 1. (color online) (a) Schematic of the device: a 
graphene ribbon is fabricated on a Si02 substrate. Source, 
drain and two side gates are contacted by gold contacts (rib¬ 
bon layer). A hBN flake is deposited on top of the device. On 
the hBN flake, an additional graphene flake is deposited and 
contacted by gold contacts (top layer), (b) False-color SFM 
image of the final device: the graphene top layer is colored 
in white and the corresponding metal contacts in red. The 
ribbon layer contacts as well as the overlay of the ribbon (not 
visible in SFM scan) are marked in blue. The orange back¬ 
ground is the hBN flake, (c) Zoom into [b]: the top layer is 
not flat on a large scale but the area directly above the ribbon 
is wrinkle-free. The wrinkles and bubbles are highlighted with 
red dashed lines, (d) SFM image of the graphene nanoribbon 
before transfer of the hBN flake. 


Figure 2. (color online) (a) Two-terminal resistance of the 
graphene top layer as a function of Vtl and Vbg- The orange 
box marks the range of the following plots, (b) Derivative of 
[a] along Vbg- The change in sign (purple line) marks the po¬ 
sition of the charge neutrality point of the top layer that is not 
above the ribbon layer, (c) Derivative of [a] along Vtl- The 
yellow line marks the position of the charge neutrality point 
in the top layer that is located above the ribbon layer. The 
black dashed lines mark the range where charge puddles are 
expected (disorder density) and the arrow marks the position 
where the kink of the line is the strongest, (d) Current flow¬ 
ing through the ribbon as a function of Vtl and Vbg voltage 
at a fixed bias voltage. The green lines mark the position of 
the region of suppressed conductance. (e,f) Absolute value of 
the calculated charge carrier density for the same range as in 
[c,d]. 


hundred Ohms are not subtracted. There are two pro¬ 
nounced features in Fig. each being roughly parallel 
to one of the gate axes. Fig. shows the derivative 
of the resistance along Vbg- At the position where the 
derivative changes sign, a slope of one (purple line) is 
observed. This line corresponds to the charge neutrality 
point of the parts of the top layer graphene that are ex¬ 
ceeding the area of the ribbon layer graphene. Therefore, 
the density is kept constant if Vtl is compensated by an 
equal amount of Vbg- In the following, the term “top 
layer” will only be used for the graphene top layer part 
lying above the ribbon layer. 


The derivative along Vtl is more complex as shown in 
Fig.^: the line along which the derivative changes sign 
(yellow line) - i.e. the position of the charge-neutrality 
point of the top layer covering the ribbon layer - is 
not straight but bent. This is expected for such a sys¬ 
tem of stacked graphene layers m and the position of 
the kink (black arrow) marks the position of the aver¬ 
age charge neutrality point of the ribbon layer. At the 
kink, the average charge carrier density in the ribbon 
layer approaches zero and screening of the back-gate is 
reduced. In addition to the position of the charge neu- 
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trality point of the top layer (yellow line), the outer 
two dashed lines mark the boundaries of the approxi¬ 
mate region where electron-hole puddles are expected in 
the graphene top layer due to disorder originating from 
fabrication. These values are extracted at a constant 
Vbg- For sufficiently high density in the ribbon layer, 
values of n^is « ±1 x 10^^ cm“^ are found for the top 
layer independently of the chosen Vbg ES] • A density of 
±1 X 10^^ cm ^ corresponds to a mean spacing of about 
30 nm between charge carriers inside the graphene layer, 
and is therefore larger than the distance between the two 
layers (13nm). 

Fig.[2ji shows the current flowing through the ribbon 
as a function of Vtl and Vbg at a DC-bias of ImV. 
Unexpectedly, there are two regions of suppressed con¬ 
ductance (blue, compare e.g. |T71[T^Eni)- The right one 
coincides well with the bending of the yellow line (see 
arrow) and can therefore be interpreted as the average 
charge neutrality point of the ribbon layer. The origin 
of the left region of suppressed conductance is unclear - 
we believe that it might originate from inhomogeneous 
doping due to fabrication residues. 

In Figs.[^,f, the average densities expected in both the 
top (nr) and ribbon layer (nji) graphene are calculated 
from [IIIET]: 


Vsa - + £,(„ri/e 


C; 


Si02 


( 1 ) 


CTI'J' 

Vtl — Vtl.D irac = EF{nii)/e — Ep{nT)/e — —— (2) 

gbn 


EF{n) = - — rhvF^/^ 


(3) 


Vtl.D irac = -0.05 V and Vbg,D irac = 7.5 V denote 
the values where both the ribbon layer and the top 
layer are simultaneously at the charge neutrality point. 
vf = 10® m/s is the Fermi velocity of charge carri¬ 
ers in graphene. The capacitances are calculated based 
on a parallel plate capacitor model (with €si 02 ~ 3.9, 
ehBN « 4, thBN « 13nm, tsi 02 « 285nm). For both 
layers, the evolution of the Dirac point in Vbg and Vtl is 
generally well predicted - despite the simplification of the 
parallel plate capacitor model which is, strictly speaking, 
not valid for the ribbon itself. Also the idealized Epin) 
relation will be modified in experiment by disorder as will 
be discussed later in more detail. 

Next, the transport regime of suppressed conductance 
in the ribbon is investigated (so-called “transport gap”). 
Figs. ii ,b show the current flowing through the ribbon 
as a function of Vbg and bias voltage for two different 
values of Vtl- In agreement with previous work, par¬ 
tially overlapping Coulomb blockade diamonds are ob¬ 
served [HHIllEIlES], indicating the presence of multiple 
sites of localized charge inside the ribbon. Figs.|^,d show 
the current at fixed bias flowing through the ribbon in 
the right transport gap as a function of the applied side 
gate voltages (fixed Vbg and Vtl)- While details vary, 
predominantly slopes of -1 are found indicating that the 
localized charges are equally well coupled to both side 
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Figure 3. (color online) (a,b) Current flowing through the rib¬ 
bon as a function of applied bias and Vbg for the right region 
of suppressed conductance. (c,d) Current flowing through 
the ribbon as a function of applied side gate voltages at a 
fixed Vbg and Vtl- For better visibility, the features of en¬ 
hanced current are extracted and shown in the middle. (a,c) 
are recorded at low and (b,d) at high top layer charge carrier 
density. 


gates and therefore likely located in the ribbon. Addi¬ 
tionally, other slopes are found indicating states that are 
not positioned symmetric relative to the ribbon axis and 
might be localized along the edges [55]. The observation 
of avoided crossings between different resonances indi¬ 
cates strong capacitive coupling between different sites 
of localized charge [29] . No significant difference is found 
when comparing transport in the right transport gap for 
different densities in the top layer (compare Figs. ii,b 
and Figs.[^,d). The behavior of the left transport gap is 
qualitatively similar to the right transport gap but cur¬ 
rent is generally less suppressed and the Coulomb block¬ 
ade diamonds are smaller. 

In the following, the evolution of the Coulomb reso¬ 
nances as a function of Vbg and Vtl is investigated as 
shown in Fig. j^. In order to enhance the visibility of 
the features, the axes were transformed (details see fig¬ 
ure caption). It is instructive to interpret the bending 
of the lines as the influence of the quantum capacitance 
rather than screening (both interpretations are equiva¬ 
lent). Eqs. 01 and 1 ^ can be rewritten using the defini¬ 
tion of the quantum capacitance Cq |27] : 

CQ{n) = e^D{EF) = ^^^^ (4) 
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Figure 4. (color online) (a) Current flowing through the ribbon as a function of Vtl and Vbg at a hxed bias. In order to improve 
the electrical stability of the measurements, several rectangular regions in Vbg and Vtl were measured consecutively. For each 
of these measurements, the axes were transformed (rotation followed by horizontal squeeze - similar to the n-D transformation 
in bilayer graphene [30]) in order to obtain a better visibility of the features. Finally, the different measurements were put 
together in one plot. Lines of constant Vtl and Vbg are marked in green and brown. The direction of the axes Vtl and 
Vbg is marked with black arrows, (b) Same data as [a] where a number of different Coulomb resonances are marked, (c) 
Slopes extracted from the lines in [b] as a function of the applied gate voltages, (d) Local quantum capacitance of the top 
layer calculated from the slopes in [c] using different thicknesses for the hBN layer, (e) Local charge carrier density of the top 
layer using the data from [d] and the same hBN thicknesses as before. Magenta curves are shifted to the left, (f) Calculated 
“capacitance disorder” ric^dis with a locally varying dielectric thickness of 2nm and a mean thickness of 14.5 nm. The brown 
lines mark the position at which this “capacitance disorder” starts to dominate over the intrinsic fabrication disorder of around 
10“cm~^ 


BVtl 

dVBG 


riji—const 


BVtl 

dVsG 


ut— const 


Csi02 _ Csi02 
Cbn 
C si02 


Csi02 + GQinpi) 


(5) 

( 6 ) 


The slope of features of constant density in the rib¬ 
bon layer {dVTh/dVBG \nR=const) yields direct informa¬ 
tion about the quantum capacitance of the top layer 
{CQ{nT))- As Coulomb resonances in the ribbon mark 
the position where one electron is added to the ribbon, 
they are used to extract such slopes as shown in Fig. [^. 


Fig.|4): shows the value of the slopes as a function of gate 
voltages. It is difficult to exactly track the resonances due 
to avoided crossings between different resonances lead to 
small wiggles in the extracted slopes. There are at least 
3 sets of slopes that are marked in different colors. Due 
to the localization of charge carriers in the ribbon, the 
Coulomb resonances are a local sensor for the electronic 
properties of the top layer. Also note that Eqs. an m 

are generic and do not contain graphene specific terms. 

We now use the evolution of the Coulomb blockade 
resonances to obtain information on the local quantum 
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capacitance of the top layer. In the limit of high top 
layer density (left side of Fig. |4:, saturation of slope), 
the quantum capacitance is negligible and the slope is 
given by —Csi 02 /CBN- We find a different saturation 
value of dVTL/dVBG for different sets of lines. While it is 
possible that different sites of localized charge experience 
a different capacitance ratio due to geometry, we believe 
it to be more likely that the observed differences originate 
from a local change in Cbn due to fabrication residues 
that locally increase the thickness of the dielectric. 

Aside from a different saturation value, the different 
curves are shifted horizontally, indicating a different lo¬ 
cal doping of the top layer graphene that is probed by 
the localized state in the ribbon. The shift between the 
three different sets of lines can be explained by different 
induced densities due to the different local capacitances 
together with an nonzero intrinsic doping of the top layer 
graphene. 

As found in previous experiments, localized states can 
slightly change their spatial position when an additional 
electron is added or when the electrostatic environment 
is changed [HIST]. The deviations within one set of lines 
can therefore be explained by slightly different areas of 
the top layer that are probed by the localized states in 
the ribbon. 

In Fig. 1^, the extracted slopes are converted into a 
local quantum capacitance of the top layer using Eq. 
and different hBN thicknesses based on the saturation 
value of the slope. The quantum capacitance is converted 
into a local charge carrier density using Eq. and by 
shifting the magenta curves by 5V to the left as shown 
in Fig. 1^ (the yellow lines are not shown for clarityh 
The thick brown line {ntheory) is calculated using Eq. 0 
together with the assumption that the middle layer is at 
zero density. 

Generally, the lines extracted from the measurement fit 
well to the expected values. The blue lines have a much 
wider spread around the theoretical value than the red 
lines. This is compatible with the model of a spatially 
changing capacitance between the ribbon and the top 
layer due to fabrication residues: the red lines correspond 
to areas that are mostly flat, whereas the blue lines corre¬ 
spond to areas with residues. As the areas on which the 
charges are localized in the ribbon are likely larger than 
the typical area of fabrication residues EH HE], spatial 
averaging of the capacitance occurs. This is compatible 
with the fact that the extracted value from the measure¬ 
ment (« 1 nm) is smaller than the value from the SFM 
scan (5nm). The shift of the magenta lines indicates a 
local difference in doping. 

In order to estimate in which regime this local change 
of capacitance is important, Eqs. m and are used to 
calculate the difference in density (jic^dis) if the thickness 
of the hBN flake is varied locally by 2nm as shown in 
Fig. Due to the plate capacitor geometry, ric^dis is 
equally large for both ribbon and top layer. As soon as 
the brown lines are crossed (uc.dis = ndis), the effect of 


varying dielectric thickness starts to dominate over the 
intrinsic doping fluctuations due to fabrication residues. 
These conclusions are important, for example, for double¬ 
gated bilayer graphene devices where a high top-gate 
voltage is applied in order to open a bandgap [52HES] . 
Choosing a thick dielectric together with as clean inter¬ 
faces as possible mitigates these layer thickness fluctua¬ 
tions. 

Finally it is possible to extract the average quan¬ 
tum capacitance of the ribbon layer using the data from 
Fig. Hi together with Eq. and tracking the position 
of the Dirac point in the top layer. The extracted 
density as a function of the gate voltages fits reason¬ 
ably well to the expectations and a disorder density 
ndis,R « I — 2 X 10^^ cm“^ is obtained for the ribbon 
layer. 

In summary, we have investigated electrical transport 
through a graphene nanoribbon and a strongly capaci- 
tively coupled graphene top layer. As this experiment is 
a crucial step towards stacked - and therefore strongly 
coupled - graphene nanostructures, it is important to un¬ 
derstand transport in this system in detail. We showed 
that the mutual interactions can be discussed in terms 
of screening or, equivalently, quantum capacitance. In 
agreement with previous studies, transport through the 
nanoribbon was governed by Coulomb blockade in the 
regime of low density. The evolution of single Coulomb 
blockade peaks was tracked as a function of voltages ap¬ 
plied to the back gate and the top layer graphene. It 
was shown that the position of the Coulomb resonances 
in gate voltage space is directly related to the quantum 
capacitance and therefore, the density of the top layer 
graphene sheet. A detailed analysis of the different ex¬ 
tracted Coulomb peak positions showed different geomet¬ 
rical capacitances per unit area between localized charges 
in the ribbon and the top layer graphene. This be¬ 
haviour can be explained with fabrication residues chang¬ 
ing the thickness of the dielectric locally. These find¬ 
ings show experimentally that it is crucial for stacked 
graphene devices to have interfaces that are as clean 
as possible. Small variations in dielectric thickness can 
quickly become more important to transport than the 
intrinsic electron-hole puddles originating from nonzero 
chemical doping of the devices. Finally, it is worth not¬ 
ing that charge localization inside the graphene nanorib¬ 
bon is not influenced strongly by the density of the top 
layer graphene, except for the mentioned evolution in 
gate voltage. 
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